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ABSTRACT 


Sound waves propagating through the oceans are refracted by internal waves. In 
the shallow waters of the continental shelves, an additional downward refraction of sound 
waves due to internal waves can cause them to interact more often with the seabed, 
resulting in additional energy from the sound waves being dissipated into the seabed. 

This study investigates how internal waves affect sound propagation on the 
continental shelves. It first quantified the types of internal waves on the continental shelf 
of California, near Point Sal, using data collected from a field experiment. Next, the 
effects that these internal waves have on sound propagation were quantified via 
simulations using a ray theory acoustic model. 

The results showed that internal waves in the experiment area were largely 
generated by tidal forcing. Compared to simulations without internal waves, simulations 
accounting for the effects of internal waves resulted in higher sound energy loss, as 
internal waves tend to cause sound waves to strike the seabed at steeper angles and over 
shorter distances. Thus, to enable a more accurate assessment of underwater acoustic 
system performance, the effects of internal waves on sound propagation in shallow 
waters need to be accounted for. 
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I. INTRODUCTION 


Internal waves have been observed in many parts of the world’s oceans (Jackson 
2004, see Figure 1), and sound waves propagating through the ocean are refracted by 
these internal waves. In deep waters, such refraction generally only causes variability in 
transmission time (Marwoto 2015); it does not cause the sound waves to have additional 
sea-surface/seabed boundary interaction given the depth of the water column. 


Figure 1. Locations of Observed Internal Waves. Source: Jackson (2004). 



The red dots indicate the locations where internal waves are recorded in “An Atlas of 
Oceanic Internal Solitary Waves.” The lack of observed internal waves in the open ocean 
is due to a lack of data at those locations. 


In shallow waters such as the continental shelves however, an additional 
downward refracting of sound waves due to internal waves could cause them to interact 
more often with the seabed, causing greater transmission loss than would be expected in a 
scenario without internal waves, as additional energy from sound waves gets dissipated 
into or gets scattered by the seabed. Apel et al. (1997, 2007) have shown the importance 
of internal waves in shallow water acoustics, while Zhou et al. (1991) have shown 25 dB 
signal fluctuations in the Yellow Sea with timescales of minutes to hours. 
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The objectives of this thesis are to identify the types of internal waves in the 
shallow waters of a continental shelf and to analyze the effects they have on sound 
propagation ranges. The internal wave data were obtained from an Office of Naval 
Research (ONR)-sponsored experiment titled “2015 Point Sal Inner-shelf Experiment 
(PSIEX),” which was conducted in the coastal waters off Point Sal, California. The 
internal waves were measured using vertical strings of temperature sensors, while their 
effects on sound propagation ranges were quantified via simulations using a ray theory 
acoustic model called BELLHOP. The end goal is to do a generic study to determine how 
internal waves affect sound propagation ranges in the shallow waters of the continental 
shelves. 
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II. PHYSICAL BACKGROUND 


A. WHAT ARE INTERNAL WAVES? 

Internal waves (IW) are gravity waves that oscillate in water, rather than on the 
sea surface. An analogy would be a layer of oil sitting on a layer of water—a small shake 
of the mixture will create waves at the oil-water interface. In order for internal waves to 
exist there must be an increase in sea-water density with increasing depth, so that any 
perturbation would then cause an underwater wave to ripple outward as internal waves. 

B. INTERNAL WAVES MANIFESTED AS WATER TEMPERATURE 
CHANGES 

In the linear equation of state, density (p) is inversely proportional to temperature 
(7) and proportional to salinity (5). 

Linear Equation of State: ——— = -a{T - T^) + P{S - S^), 

Po 

where po =1025 kg/m^, So = 35p su, To=\0 °C, a = 2x10'"^ °C'', P = 7.5x10“^; a is the 
thermal expansion coefficient and P is the coefficient of saline contraction. 

For this study, a constant salinity throughout the water column was assumed, as 
there was no introduction of fresh water into the region. Hence, temperature changes 
were the only factor contributing to changes in density in this investigation of internal 
waves. As deeper waters are usually colder than shallower waters, this implies that they 
are also denser than shallower waters. Thus, any external forcing would set up a wave at 
the interface between the waters of different temperature (and density). In other words, an 
internal wave, which propagates as perturbations in the density difference of water, is 
manifested and can be recorded as changes in water temperature. 
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c. 


TYPES OF INTERNAL WAVES IN THE SHALLOW WATERS OF THE 
CONTINENTAL SHELVES 


1. Internal Tides—Internal Waves Generated by Tidal Forcing 

Gerkema and Zimmerman (2008) show that internal waves on the continental 
shelves are generated by tidal forcing over the depth transition between the deep ocean 
and the shallow continental shelves. In this scenario, the internal waves have a semi¬ 
diurnal frequency of two cycles per day (cpd). This agrees well with our experiment 
observations where over a period of one day, two internal waves due to tidal forcing were 
observed (Figure 2). Internal waves generated by tidal forcing are also called internal 
tides. 

Very often, a steep interface can occur between the cold and warm front of an 
internal wave. Such a steep interface is called a tidal bore or, simply, bore (see Figure 2). 


Figure 2. Types of Internal Waves Observed in One Day. 
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2. Solitons—Solitary Internal Waves 

In addition to internal tides, shorter duration internal waves were observed 
throughout the data. These short spikes of internal waves are called internal solitary 
waves or solitons (Figure 2), which are generated from internal tides transferring their 
energies to shorter-space scales (Gerkema and Zimmerman 2008). 
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D. SEABEDS AND THEIR EEFECTS ON SOUND PROPAGATION 


The types of seabed in shallow waters play a primary role in determining the 
range of sound propagation, as the waters in these areas are usually downward refracting 
and will cause sound to interact with the seabed multiple times over short distances 
(Jensen et al. 1994). Figure 3 illustrates this effect. 


Figure 3. Ray Trace of a Horizontally Launched Ray in Mid-water Column. 



Left: A typical slightly downward refracting sound speed profile recorded during the 
experiment was used in the ray trace simulation run as shown to the right. 


Right: A mid-water column horizontally launched sound “ray” gets refracted downward 
and interacts with the seabed eight times in 10 km. 


Three representative continental shelf seabed types were used in the simulation 
runs: rocky seabed (chalk), sandy seabed (sandy-slit), and muddy seabed (clay-mud). The 
detailed parameters of each bottom characteristic used in the BELLHOP runs are listed in 
Table 1. 


Table 1. Bottom Types Parameters Used in BELLHOP Runs. 

Adapted from Colosi (2016). 


Bottom Type 

Compressional 
sound speed 
(m/s) 

Density 

(kg/m') 

Attenuation 

(dB/X) 

Rocky (chalk) 

2400 

2200 

0.2 

Sandy (sandy-silt) 

1644 

1769 

1.22 

Muddy (clay-mud) 

1470 

1312 

0.09 
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To understand the effeets of bottom interaction on acoustics energy, it is 
necessary to invoke the formula for reflection coefficient to determine how much of the 
original ray intensity was reflected from the seabed at various grazing angles. 


m sin 9 -Jn^- cos^ 9 

Reflection coefficient: = - -—, 

m sin 9^ + Jn^ - cos^ 9^ 

where 9g is the sound ray grazing angle as measured from the seabed plane, 
P C 

m = — , n = —, p\ and p 2 are seawater and seabed density, respectively, and c\ and C 2 
P\ ^2 

are sound speed in seawater and seabed, respectively. 

It should be noted that while all the seabed types are denser than seawater 
(approximately 1025 kg/m ), the sound speeds in rocky and sandy seabed are faster than 
seawater (approximately 1500 m/s), while the sound speed of a muddy seabed is slower 
than that of seawater. This has important physical implications—for the rocky and sandy 
seabed, total internal reflection (|Ri 2 | = 1) occurs when the direction of sound propagation 
(or a sound “ray”) strikes the seabed at angles, 9g < 9cr, where 9cr (critical angle) is the 
highest angle at which total internal reflection can occur (Figure 4). As for the case of a 
muddy seabed, total internal reflection does not occur as there is no critical angle. Instead 
there is what is called a Brewster or intromission angle where |Ri 2 | = 0, in which all 
energy from the sound waves in water gets dissipated into the seabed when 9g is at the 
Brewster angle. 
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Figure 4. Reflection Coefficients of Rocky, Sandy, and Muddy Seabeds for 

the Various Grazing Angles. 



Bottom Type 

— Rocky (chalk) 

— Sandy (sandy-silt) 
— Muddy (clay-mud) 



A sound “ray” is a vector that is normal to the sound wave front and which traces the path 
that a sound wave travels. 


Referring to Figure 4, the rocky seabed has the highest critical angle of 
approximately 52°; hence, any sound ray striking a rocky seabed at less than that angle 
will experience total internal reflection. In fact, more than 55% of the sound energy will 
be reflected even when sound strikes the seabed at a 90° grazing angle. This means that 
the rocky seabed will allow sound to travel the farthest as it is the least lossy of the three 
seabeds. 

For a sandy seabed, the critical angle is much smaller compared to the rocky 
seabed at approximately 25°. At grazing angles greater than the critical angle, the portion 
of reflected energy quickly drops to below 40% at grazing angles of above 40°. Hence, a 
sandy seabed is more lossy compared to a rocky seabed, and sound experiences more 
attenuation when it interacts with a sandy seabed. 

For the muddy seabed, the percentage of energy reflected drops very rapidly from 
100% at 0° grazing angle to 20% at 6° grazing angle. The Brewster angle is reached at 
approximately 12.5° grazing angle, and all the energy from the sound waves gets 
absorbed into the seabed. Beyond the Brewster angle, up to 90° grazing angle, only 
approximately 12% of the sound energy gets reflected back into the water. Hence, a 
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muddy bottom is extremely lossy, and sound will experience the most attenuation in this 
environment. 

The simulated transmission loss plots for the three seabed types are as shown in 
Figure 5. The plots were generated using the sound speed profile (SSP) of Figure 3, and 
the only difference between the three plots is the seabed type. From the transmission loss 
plots, it is apparent that a rocky bottom has the least sound energy transmission loss, 
while a sandy bottom has more transmission loss and a muddy bottom has the most 
transmission loss. 

Figure 5. Transmission Loss Plots for a Rocky Seabed (top), 

a Sandy Seabed (center), and a Muddy Seabed (bottom). 
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III. EXPERIMENT OVERVIEW 


The data used in this investigation was collected during an ONR-sponsored effort 
titled “2015 Point Sal Inner-shelf Experiment (PSIEX).” The objective of the experiment 
was to quantify the various oceanographic phenomena in the coastal inner-shelf. The 
experiment took place in the waters off Point Sal, California, during the summer of 2015. 

The locations of the various sensors are as shown in Figure 6. The experiment 
lasted 52 days from 15-Jun-2015 07:00 PDT (yearday 166.29) to 06-Aug-2015 07:00 
PDT (yearday 218.29). 


Figure 6. Experiment Site Sensor Deployment Locations. 

Source: Colosi et al. (2016). 
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For this study, the focus was on temperature string 050, as it was the deepest 
deployed sensor of the experiment at 52 m water depth. Temperature string 050 
consisted of five temperature sensors deployed from a depth of 6.2 m to 46.1 m at 
approximately 10 m spacing, and the water temperature was sampled once every 30 
seconds (or 120 samples per hour). This allowed an internal wave frequency of up to 60 
cph to be sampled, based on Nyquist-Shannon sampling theorem, which states that for a 
waveform to be accurately sampled, the sampling frequency must be two times that of the 
sampled waveform’s highest frequency. 

The maximum vertical oscillating frequency of internal waves is limited by the 
Brunt-Vaisala frequency, or buoyancy frequency, and any external forcing of a water 
column at a frequency that is higher than the Brunt-Vaisala frequency will not generate 
internal waves. Robert and William (1983) observe that the typical Brunt-Vaisala 
frequency in shallow water is 12-15 cph. Thus, the temperature sensor sampling rate of 
120 samples per hour is adequate for observing the internal waves expected in this 
environment. 
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IV. CHARACTERISTICS OF INTERNAL WAVES IN THE 

EXPERIMENT AREA 


A. OVERVIEW OE THE ENTIRE TEMPERATURE TIME SERIES 

The recordings of temperature string 050 as a function of depth and time are 
shown in Figures 7a-h, with Figures 7a-g each showing one week of data and Figure 7h 
showing three days of data from week 8. 


Figure 7. Temperature Plots of Temperature String 050 for the 
Entire Experiment (Figures 7a-h). 


Figure 7a. Temperature Plot—Week 1. 
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Figure 7b. Temperature Plot—Week 2. 
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Figure 7c. Temperature Plot—Week 3. 
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Figure 7d. Temperature Plot—Week 4. 
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Figure 7e. Temperature Plot—Week 5. 
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Figure 7f. Temperature Plot—Week 6. 
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Figure 7g. Temperature Plot—Week 7. 
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Figure 7h. Temperature Plot—Week 8, First 3 Days. 
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B. INTERNAL TIDES 

From these plots, it was observed that there was great variability in internal wave 
activity, ranging from internal waves with small temperature difference (e.g., yeardays 
174-177), to internal waves with large temperature difference (e.g., yeardays 204-208). 
The internal waves were observed to have two oscillations per day, matching the 
semidiurnal tidal frequency. Hence, it can be concluded that internal waves in this region 
were largely due to tidal forcing. 

C. SOLITONS 

Figure 8 shows the number of soliton occurrences per day. Great variability in the 
daily occurrence of solitons was also observed, ranging from no solitons during yeardays 
174 and 178 to 37 solitons during yearday 210. 
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Figure 8. Daily Occurrence Rate of Solitons. Source: Colosi et al. (2016). 
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D. INTERNAL WAVE PERIOD, PROPAGATION SPEED, AND DIRECTION 

Cross-correlation of temperature-strings in the experiment area (Colosi et al. 
2016) showed that the internal waves propagated west to east at an average speed of 0.25 
m/s (Figure 9, top), and at an average bearing of 80° with respect to true north (Figure 9, 
bottom). For ease of illustration, a purely eastward propagation of internal waves would 
be assumed for the study, as this is a general analysis of how internal waves affect sound 
propagation and is not tied to a specific bathymetry. 
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Figure 9. Speed of Internal Wave Propagation (top) and Direction of Internal 
Wave Propagation (bottom). Source: Colosi et al. (2016). 



Internal tides in shallow waters are also called tidal bores, or simply bores. The angles 
used in the “Bore Direction” plot are with respect to east, and positive angles are anti¬ 
clockwise from east. 


E. INTERNAL WAVE FRONTS IN THE EXPERIMENT AREA 

The internal waves will be assumed to have a plane wave structure as they 
propagate eastward. This implies that sound waves that travel eastward or westward will 
cut across the internal wave fronts (red arrow of Figure 10); hence, a range-dependent 
SSP was used to simulate the behaviors of these sound waves. As for the sound waves 
that travel northward or southward, they will travel in-parallel or along the internal wave 
fronts (green arrow of Figure 10); hence, a range-independent sound speed profde was 
used to simulate the propagation behavior of these sound waves. 
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Figure 10. Sound Propagation Along and Across Internal Wave Fronts in the 

Experiment Area (figure not drawn to scale). 
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Assuming that the internal waves propagate eastward and have a plane wave structure, 
easterly/westerly sound propagation will cut across internal wave fronts (red arrow) and 
experience range-dependent sound propagation, while northerly/southerly sound 
propagation will run along internal wave fronts (green arrow) and experience range- 
independent sound propagation 


F. TRANSLATING WATER TEMPERATURE TO SOUND SPEED 

For the translation of temperature, depth, and salinity into sound speed, the nine- 
term Mackenzie equation was used (Mackenzie 1981), where sound speed (c) is a 
function of temperature (T) in degrees Celsius, salinity (S) in psu, and depth (D) in 
meters. 
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For the nine-term Mackenzie equation: 

c(A‘S',7’) =1448.96 + 4.5917’ - 5.304 x 10""7’" + 2.374 x lO-^r'+ 1.340 {S-35) + 
1.630 X 10'"79 + 1.675 X 10"’79"- 1.025 X 10-"r(5' - 35) - 7.139 x lO'*"779" 

the range of the equation validity is temperature: 2 to 30°C; salinity: 25 to 40 psu; depth: 

0 to 8000 m. 

From the Mackenzie equation, it can be seen that in shallow waters, changes in 
temperature are the primary cause of sound speed changes, while sound speed changes 
due to salinity changes are small in the experiment area as there was no introduction of 
fresh water into the region. An average salinity of 34.44 psu was used for the simulation 
runs. 

Translation of wave measurements in time into distance traveled could be 
obtained by multiplying measurement time with internal wave propagation speed of 0.25 
m/s. For semidiurnal internal tides (approximately 12 hours apart crest-to-crest), this 
translates to about 10.8 km between the crests of the internal tides. Over two days, the 
internal waves will propagate 43.2 km. 

With the translation of internal wave measurements in time to physical distance, 
and the Mackenzie equation to get the SSP, there were sufficient data to generate a range- 
dependent sound speed profile to analyze the effects of internal waves on sound 
propagation. 

An example of translating temperature readings into a range-dependent SSP is 
shown in Figure 11 for the two-day period of yeardays 200-202. The temperature as a 
function of time and depth plot of this particular period is as shown in the top plot of 
Figure 11. From this data, a range-dependent SSP (with depth adjustments to facilitate 
BELLHOP analysis) was derived, as shown in the bottom of Figure 11. For the x-axis, 
the time measurements were converted to range based on the internal wave propagation 
speed, with 0 km at yearday 200, and yearday 202 converts to 43.2 km west of 
temperature string 050 (as the internal wave fronts at that location will reach temperature 
string 050 two days later with an eastward propagation). Such range-dependent SSPs 
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were fed into BELLHOP for simulations runs, to determine their effects on sound 
propagation. 


Figure 11. Temperature as a Function of Time and Depth (top) 

Translated to a Range-Dependent Sound Speed Profde (bottom). 
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It should be noted that an increase in temperature leads to an increase in sound 
speed; hence, the structure of the temperature plot and the range-dependent SSP plot look 
similar. This means that the structure of the internal waves has a direct influence on 
sound speed and the shape of the internal waves can be seen from the range-dependent 
SSP. 
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V. BELLHOP—THE MODELING AND SIMULATION MODEL 

CHOSEN 


A. RAY TRACING 

A ray theory acoustic model called BELLHOP was used to ray trace the sound 
wave propagation path through the water. An example of a ray trace plot has been 
previously illustrated in Figure 3. The physics behind ray tracing is Snell’s Law, which 
calculates the refraction of sound waves as they move horizontally through waters of 
different density and hence sound speeds. A depiction of Snell’s Law is shown in Figure 
12 . 


Figure 12. Depiction of Snell’s Law. Adapted from Kapolka (2015). 



dcos 6 >i =A^ 
dcos6>2 =/l2 

cos 6 ) /Ij /'/ Cj 


cos (9, A, c^/ c 

' f 


c c 

—*— - —-— = constant 


cos 6 >, cos 6>2 


Cl, C 2 : Sound speed of layer 1 and layer 2, respectively. 

Wavelength of sound wave in layer 1 and layer 2. 

01, 02: Grazing angle, angle of sound “ray” with respect to the plane dividing the two 
mediums. 


c c 

Given that—^— = —-—, sound waves will refract downward as sound speed 
cos 6^ cos 02 

decreases with depth, and will refract upward as sound speed increases with depth (see 
Figure 13). 
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Figure 13. Refraction of Sound Waves. Adapted from Kapolka (2015). 


Sound 



1 2 4 


Sound 



Sound speed increasing with depth 



Top: Sound wave refracted downward as sound speed decreases with depth. 
Bottom: Sound wave refracted upward as sound speed increases with depth. 
Cl, C 2 , C 3 , C 4 : Sound speed of the various layers. 


To generate the transmission loss plots as seen in Figure 5, BELLHOP ray traces 
a large number of sound “rays” emanating from a sound source and then sums the 
intensity of the rays as they move downrange. 

B. WHY BELLHOP RAY THEORY ACOUSTIC MODEL 

A ray theory acoustic model was chosen as this study uses ray tracing of the 
sound waves to demonstrate how internal waves modulate their interactions with the 
seabed. The other reason was that a ray theory acoustic model could also support range- 
dependent SSPs (Etter 2013), which are required in this study as internal waves will 
cause the SSPs to vary with range. 
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C. BELLHOP’S ENVIRONMENTAL PARAMETERS 

1. Adjustment of Actual Sea Depth to BELLHOP Model’s Depth 

As there were no temperature sensors at the sea surface and seabed for 
temperature string 050, a translational step was made to ensure compatibility with 
BELLHOP by treating the shallowest sensor as being at the top boundary of the water 
column and the deepest sensor as being at the bottom boundary of the water column in 
BELLHOP. This eliminates the need to extrapolate the temperature readings to the sea 
surface and the seabed. 

The actual depth of the temperature sensors and their adjusted depth in 
BELLHOP are as listed in Table 2. The depth adjustment was done by treating the 
shallowest sensor depth (at 6.2 m) as the “surface sensor’’ (of 0 m) and adjusting all the 
sensor depths upward by 6.2 m, with the depth of the deepest sensor rounded to 40 m. 


Table 2. Temperature String 050—Temperature Sensor Depths. 


Actual sensor depth 

Adjusted depth for 
BELLHOP run (by -6.2 
m) 

6.2 m 

0.0 m 

16.3 m 

10.1 m 

29.7 m 

23.5 m 

36.1 m 

29.9 m 

46.1 m 

Rounded to 40.0 m 


2. Flat Bathymetry and Smooth Sea Surface 

To enable a general study of the effect of internal waves on sound propagation, a 
flat seabed and a smooth sea surface were used in the simulation runs, so as not to limit 
the applicability of the analysis to a specific location. 
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3. 


Sound Source Depth 


The sound source was placed at a depth of 38 m for the simulation runs so that 
sound rays launched at small angles with respect to the horizontal plane would not be 
refracted downward too much before striking the seabed. This allows them to graze the 
seabed at a small angle, so that most of their energy will be reflected, thus enabling them 
to propagate farther, as is the case for many underwater applications. 

D. MEASURES TO MAXIMIZE SIMULATION ACCURACY 

1. Ray Trace Using a Frequency of 2 kHz 

Etter (2013) states that the geometrical acoustics approximation made for the ray 
theory model, which assumes that sound speed does not change much over one wave 
length horizontally, limits the ray acoustics to higher frequency applications (as lower 
frequency sound waves have longer wave lengths). 

Etter (2013) also gives a guideline for defining the frequency above which a ray 
acoustics model becomes suitable for use, in the following equation: 


where / is the frequency, H is the water depth, and c is the speed of sound. With a 
BELLHOP simulation depth of 40 m, and using an average sound speed of 1500 m/s, the 
formula resolves to f > 375 Hz. A frequency of 2 kHz, which is around the frequency 
range of a mid-frequency sonar, was used for the study and is well above the 375 Hz 
threshold. 


2. Gaussian Beam Option 

Ray trace models are unable to approximate the effects of caustics, which are 
regions of intense beam focusing. One method to mitigate this is Gaussian beam tracing. 
Porter and Bucker (1987) state that Gaussian beam tracing avoids the traditional ray trace 
artefact of infinite energy levels at caustics, and is attractive for high-frequency, range- 
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dependent applications. Hence, the Gaussian beam option in BELLHOP was used to 
increase the ray trace modeling accuracy. 

3. Cubic Spline Fitting of the SSP Plot 

The depth interval between sensors on temperature string 050 was 10 meters. 
Such a coarse depth resolution could cause inaccuracies in ray tracing if straight lines 
(linear interpolation) were used to join the data points. To overcome this, the cubic spline 
interpolation of the SSP was used in the BELLHOP runs, which replaced the linear 
approximation with a more realistic curved line for the SSP (see Figure 14). 


Figure 14. Comparison between Straight Line Linear Approximation and 
Cubic Spline Interpolation of the Sound Speed Profde. 



4. Incoherent Beam Option 

The BELLHOP simulation runs were conducted using the “incoherent” beam 
option, which allows a smoothed transmission loss plot of the sound energy over range. A 
“coherent” beam option would have accounted for the phase of the sound waves and 
shown the constructive and destructive interferences of the signal over range as 
distinctive peaks and null in the transmission loss plot. The difference in transmission 
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loss between an incoherent and a coherent beam option is illustrated in Figure 15. The 
incoherent beam option was chosen as this study is focused on the effects of internal 
waves on sound propagation over longer ranges, and does not require an in-depth analysis 
of the effects of constructive and destructive interferences of the signal over short 
distances of a few wave lengths. 


Figure 15. Incoherent Transmission Loss (top), 
Coherent Transmission Loss (bottom). 
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VI. TYPES OF ANALYSIS AND THEIR CONSIDERATIONS 

A. TIME PERIODS CHOSEN FOR ANALYSIS 

Five time periods with distinct internal wave characteristics were chosen from the 
experiment data for the BELLHOP simulation runs (Table 3). The first was the period at 
yeardays 175-177, which was a period of well-defined internal tides with few solitons. 
The difference between the surface and bottom water temperature in this period was 
small at 2.5°C. 

The second period was at yeardays 205-207, which was also a period of well- 
defined internal tides, but with many solitons. The surface and bottom water temperature 
difference in this period was large at 7°C. 

The third period was at yeardays 210-212, which was a period of weak internal 
tides but with a large number of soliton occurrences and a surface and bottom water 
temperature difference of 6°C. 

The fourth period, at yeardays 200.6-202.9, had well-defined internal tides but 
few solitons and a surface and bottom water temperature difference of 5°C. Two 
simulation runs were conducted for period 4, one with the sound source in the middle of 
the cold front of an internal tide (yearday 200.6, period 4a of Table 3), the other with the 
sound source in the middle of the warm front of an internal tide (yearday 200.9, period 4b 
of Table 3). 

The fifth period was at approximately yearday 174.22, which was at the warm and 
cold interface of a steep tidal bore. Two simulation runs were conducted at data points 
that were 5 minutes or 75 m apart, with one run at the warm front of the tidal bore and the 
other run at the cold front of the tidal bore. 
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Table 3. Internal Wave Charaeteristies of the Seleeted Time Periods for Analysis. 
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B. BELLHOP PLOTS USED FOR ANALYSIS 


Transmission loss plots are generated using range-independent and range- 
dependent SSPs, so that a quantitative comparison could be made on how much internal 
waves affect sound propagation (range-dependent SSP) as compared to a model that does 
not account for internal waves (range-independent SSP). An example of a range- 
independent and range-dependent sound speed profde is shown in Figure 16. 


Figure 16. Range-Independent SSP (top), Range-Dependent SSP (bottom). 
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It has been mentioned that a range-independent SSP would be more appropriate 
for simulating along internal wave front sound propagation (northward/southward in the 
experiment area), while the range-dependent SSP would be more appropriate for 
simulating across internal wave front sound propagation (eastward/westward in the 
experiment area). 

Two types of transmission plots are generated; the first is the “traditional” 
transmission loss plot that has been shown in Figure 5, which uses a gradual color 
scheme to show the gradual loss of signal strength over range. Such plots allow the sound 
loss “structure” to be seen. 

The second type of transmission loss plot is the 70 dB transmission range plot, an 
example of which is shown in Figure 17. Such plots allow a quantitative assessment of 
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the range a sound signal can reach before its transmission loss exceeds 70 dB. Referring 
to Figure 18, yellow shows the region where transmission loss is less than 70 dB, while 
dark blue shows the region where transmission loss is more than 70 dB. A transmission 
loss of 70 dB was chosen as the cut-off value in this study as it allows the display area of 
the plot to be maximized, which for the case of a rocky bottom of Figure 18, the 70 dB 
transmission loss is at 40 km, which is almost to the end of the plot. For this study, the 70 
dB transmission ranges are measured from the source range (0 km) to the maximum 
range that the 70 dB contour reaches. 


Figure 17. Example of a 70 dB Transmission Range Plot. 
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VII. RESULTS OF MODELING AND SIMULATION RUNS 


A. PERIOD 1: INTERNAL TIDES WITH SMALL TEMPERATURE 

DIFFERENCE, FEW SOLITONS 

The first period ehosen for analysis (yeardays 175-177) was eharacterized by 
well-defined internal tides, with a small temperature differenee of 2.5°C (from 10 to 
12.5°C) between the bottom and the surfaee of the water, and few solitons (fewer than 5 
solitons a day, see Figure 8). 

The sound speed profile ranged from approximately 1489 m/s to 1498 m/s. The 
range-independent sound speed profile and range-dependent sound speed profile are 
shown in Figure 18, and the SSP at 0 km is at yearday 175. It can be seen from the range- 
independent sound speed profile (Figure 19, top) that a bottom layer exists from 30-40 m 
depth, and this should serve to give an extended transmission range in the bottom layer 
due to downward refraction trapping sound energy in the lower water column. 

The transmission loss plots for the three types of seabed, with a range- 
independent SSP and range-dependent SSP, are shown in Figure 19. The 70 dB 
transmission range plots are shown in Figure 20. The sound source depth for all these 
plots is 38 m. 

Across all bottom types, the modeling runs that accounted for internal waves gave 
a reduced transmission range compared to a modeling run that does not factor the effects 
of internal waves. For the rocky seabed, the 70 dB transmission range dropped from 40.0 
km for along internal wave front propagation to 21.5 km for across internal wave front 
propagation. For the sandy seabed, the 70 dB transmission range dropped from 27.3 km 
to 6.2 km, while for the muddy seabed, the 70 dB transmission range dropped from 26.5 
km to 2.5 km. Hence, it can be seen that for a small temperature difference between the 
seabed and sea surface, and with a bottom layer, factoring the effects of internal waves 
resulted in significant acoustics range reduction compared to a model that does not 
account for internal waves. 
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Figure 18. Period 1—Range-Independent SSP (top), Range-Dependent 

SSP (bottom) 0 km at Yearday 175. 
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Figure 19. Period 1—70 dB Transmission Range for a 

Rocky, Sandy, and Muddy Seabed. (Figures 19a-c). 


Figure 19a. Period 1—Rocky Bottom Transmission Loss through Range-Independent 

SSP and Range-Dependent SSP. 
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Figure 19b. Period 1—Sandy Bottom Transmission Loss through Range-Independent 

SSP and Range-Dependent SSP. 
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Figure 19c. Period 1—Muddy Bottom Transmission Loss through Range-Independent 

SSP and Range-Dependent SSP. 
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Figure 20. Period 1—70 dB Transmission Range for a Roeky, Sandy, and 

Muddy Seabed. (Figures 20a-c). 


Figure 20a. Period 1—Roeky Bottom 70 dB Transmission Range for a Range- 
Independent SSP and a Range-Dependent SSP. 
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Period 1, rocky bottom 70dB transmission range, range-dependent SSP (across IW front) 
(Yellow: Tx loss < 70 dB; Blue: Tx loss > 70 dB) 




Figure 20b. Period 1—Muddy Bottom 70 dB Transmission Range for a Range- 

Independent SSP and a Range-Dependent SSP. 
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Figure 20c. Period 1—Muddy Bottom 70 dB Transmission Range for a Range- 
Independent SSP and a Range-Dependent SSP. 

Period 1, muddy bottom TCldB transmission range, range-independent SSP (along IW front) 
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Period 1, muddy bottom 70dB transmission range, range-dependent SSP (across IW front) 
(Yellow: Tx loss < 70 dB; Blue: Tx loss > 70 dB) 




B. PERIOD 2: INTERNAL TIDES WITH LARGE TEMPERATURE 

DIFFERENCE, MANY SOLITONS 

This period (yeardays 205-207) was characterized by well-defined internal tides, 
with a large temperature difference of 7°C (from 11 to 18°C) between the bottom and the 
surface of the water, and many solitons (more than 15 solitons per day, see Figure 8). The 
sound speed profde ranged from approximately 1494 m/s to 1515 m/s. The range- 
independent sound speed profde and range-dependent sound speed profde are shown in 
Figure 21, and the SSP at 0 km is centered at yearday 205. The transmission loss plots for 
the three types of seabed are shown in Figure 22, while the 70 dB transmission range 
plots are shown in Figure 23. 

Similar to period 1, factoring for internal waves in the model causes more 
attenuation compared to a modeling run without internal waves. For the rocky seabed, the 
70 dB transmission range dropped from 37.5 km to 18.3 km when the effects of internal 
waves were accounted for. For the sandy seabed, the 70 dB transmission range dropped 
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from 12.2 km to 6.5 km, while for the muddy seabed, the 70 dB transmission range 
actually increased from 1.8 km to 2.5 km. 


Figure 21. Period 2—Range-Independent SSP (top), 

Range-Dependent SSP (bottom) 0 km at Yearday 205. 
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Figure 22. Period 2—Transmission Loss for a Rocky, Sandy, and Muddy 

Seabed. (Figures 22a-c). 


Figure 22a. Period 2—Rocky Bottom Transmission Loss through Range- 

Independent SSP and Range-Dependent SSP. 
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Figure 22b. Period 2—Sandy Bottom Transmission Loss through 

Range-Independent SSP and Range-Dependent SSP. 
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Figure 22c. Period 2—^Muddy Bottom Transmission Loss through 

Range-Independent SSP and Range-Dependent SSP. 
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Figure 23. Period 2—70 dB Transmission Range for a Roeky, Sandy, and 

Muddy Seabed. (Figures 23a-c). 

Figure 23 a. Period 2—Rocky Bottom 70 dB Transmission Range for a 
Range-Independent SSP and a Range-Dependent SSP. 
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(Yellow: Tx loss < 70 dB; Blue: Tx loss > 70 dB) 




Figure 23b. Period 2—Sandy Bottom 70 dB Transmission Range for a 

Range-Independent SSP and a Range-Dependent SSP. 
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Figure 23c. Period 2—Muddy Bottom 70 dB Transmission Range for a 

Range-Independent SSP and a Range-Dependent SSP. 
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C. PERIOD 3: WEAK INTERNAL TIDES, MANY SOLITONS 

This period (yeardays 210-212) was characterized by weak internal tides and 
many solitons (more than 15 solitons per day, see Figure 8). The water temperature 
difference was strong at 6°C (from 11.5 to 17.5°C) between the bottom and the surface of 
the water. The sound speed profile ranged from approximately 1495 m/s to 1514 m/s. The 
range-independent sound speed profile and range-dependent sound speed profile are 
shown in Figure 24, and the SSP at 0 km is at yearday 210. The transmission loss plots 
for the three types of seabed are shown in Figure 25, while the 70 dB transmission range 
plots are shown in Figure 26. 

The simulation showed that even with internal waves only in the form of solitons, 
they still caused a range reduction as compared to a simulation that does not account for 
internal waves. For a rocky bottom, the 70 dB transmission range was reduced from 36.8 
km for along internal wave front propagation to 18.5 km for across internal wave front 
propagation. For a sandy bottom, the range was reduced from 17.2 km to 9.5 km. For a 
muddy bottom, the range was reduced from 7.5 km to 3.7 km. 
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Figure 24. Period 3—Range-Independent SSP (top), 

Range-Dependent SSP (bottom) 0 km at Yearday 210. 
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Figure 25. Period 3—Transmission Loss for a Rocky, Sandy, and Muddy 

Seabed. (Figures 25a-c). 


Figure 25a. Period 3—Rocky Bottom Transmission Loss through Range-Independent 

SSP and Range-Dependent SSP. 
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Figure 25b. Period 3—Sandy Bottom Transmission Loss through Range-Independent 

SSP and Range-Dependent SSP. 


Period 3, sandy bottom transmission loss through a range-independent SSP ir^nsmi^^ion Lo&s 



Range (km) 


Period 3, sandy bottom transmission loss through a range-dependent SSP 



0 5 10 15 20 25 30 35 40 43,2 


Range (km) 


Figure 25c. Period 3—Muddy Bottom Transmission Loss through Range-Independent 

SSP and Range-Dependent SSP. 
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Figure 26. Period 3—70 dB Transmission Range for a Roeky, Sandy, and 

Muddy Seabed. (Figures 26a-e). 


Figure 26a. Period 3—Roeky Bottom 70 dB Transmission Range for a Range- 
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Figure 26b. Period 3—Sandy Bottom 70 dB Transmission Range for a Range- 
Independent SSP and a Range-Dependent SSP. 

Period 3, sandy bottom TOdB transmission range, range-independent SSP (along IW front) 
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Figure 26c. Period 3—Muddy Bottom 70 dB Transmission Range for a 

Range-Independent SSP and a Range-Dependent SSP. 
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D. PERIOD 4: INTERNAL TIDES WITH LARGE TEMPERATURE 

DIFFERENCE, FEW SOLITONS 

Period 4 (yeardays 200.6-202.9) was characterized by well-defined internal tides, 
with a relatively large temperature difference of 5.5°C (from 10.5 to 16°C) between the 
bottom and the surface of the water, and few solitons (fewer than 5 solitons a day, refer to 
Figure 8). The sound speed profile ranged from approximately 1492 m/s to 1508 m/s. For 
period 4a (yeardays 200.6-202.6), the sound source was placed at the center of a cold 
front of an internal tide (Figure 27), while period 4b (yeardays 200.9-202.9) had the 
sound source placed in the center of the warm front of the same internal tide (Figure 30). 
The objective was to determine if there are significant differences in sound attenuation 
when a sound source is in the different temperature fronts of an internal tide. 

For period 4a, the transmission loss plots for the three types of seabed are shown 
in Figure 28. The 70 dB transmission range plots are shown in Figure 29. For a rocky 
bottom, the 70 dB transmission range was reduced from 39.2 km for along internal wave 
front propagation to 18.3 km for across internal wave front propagation. For a sandy 
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bottom, the range was reduced from 19.5 km to 6.7 km. For a muddy bottom, the range 
was reduced from 4.3 km to 2.0 km. 

For period 4b, the transmission loss plots for the three types of seabed are shown 
in Figure 31. The 70 dB transmission range plots are shown in Figure 32. For a rocky 
bottom, the 70 dB transmission range was reduced from 37.8 km for along internal wave 
front propagation to 19.2 km for across internal wave front propagation. For a sandy 
bottom, the range was reduced from 10.0 km to 5.8 km. For a muddy bottom, the range 
remains the same at 2.0 km. 

It can be observed that for a sandy and a muddy seabed, having a sound source in 
the cold or warm front of an internal tide results in a large variability in along internal 
wave front sound propagation ranges. For a sandy bottom, the along cold front 70 dB 
transmission range was 19.5 km, while the along warm front 70 dB transmission range 
was 10.0 km. For a muddy bottom, the along cold front 70 dB transmission range was 4.3 
km, while the along warm front 70 dB transmission range was 2.0 km. 


Figure 27. Period 4a—Range-Independent SSP (top), Range-Dependent SSP 
(bottom) 0 km at Yearday 200.6 (Cold Front of IW). 
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Figure 28. Period 4a—Transmission Loss for a Rocky, Sandy, and Muddy 

Seabed. (Figures 28a-c). 


Figure 28a. Period 4a—Rocky Bottom Transmission Loss through Range-Independent 

SSP and Range-Dependent SSP. 
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Figure 28b. Period 4a—Sandy Bottom Transmission Loss through Range-Independent 

SSP and Range-Dependent SSP. 
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Figure 28c. Period 4a—Muddy Bottom Transmission Loss through Range- 

Independent SSP and Range-Dependent SSP. 
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Figure 29. Period 4a—70 dB Transmission Loss Range for a Rocky, Sandy, and 

Muddy Seabed. (Figures 29a-c). 


Figure 29a. Period 4a—Rocky Bottom 70 dB Transmission Range for a Range- 
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Figure 29b. Period 4a—Sandy Bottom 70 dB Transmission Range for a Range- 

Independent SSP and a Range-Dependent SSP. 
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Figure 29c. Period 4a—Muddy Bottom 70 dB Transmission Range for a Range- 

Independent SSP and a Range-Dependent SSP. 
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Figure 30. Period 4b—Range-Independent SSP (top), Range-Dependent SSP 
(bottom) 0 km at Yearday 200.9 (Warm Front of IW). 
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Figure 31. Period 4b—Transmission Loss for a Rocky, Sandy, and Muddy 

Seabed. (Figures 31a-c). 


Figure 31a. Period 4b—Rocky Bottom Transmission Loss through Range-Independent 

SSP and Range-Dependent SSP. 
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Figure 3 lb. Period 4b—Sandy Bottom Transmission Loss through Range-Independent 

SSP and Range-Dependent SSP. 
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Figure 31c. Period 4b—Muddy Bottom Transmission Loss through Range- 

Independent SSP and Range-Dependent SSP. 


Period 4b, muddy bottom transmission loss through a range-independent SSP Transmission Loss 



0 5 10 15 20 25 30 35 40 43.2 

Range (km) 

Period 4b, muddy bottom transmission loss through a range-dependent SSP 
Across internal wave front sound propagation 



0 5 10 15 20 25 30 35 40 43.2 

Range (km) 


47 












Depth(m) Depth(m) Depth{m) Depth(m) 


Figure 32. Period 4b—Transmission Loss Range for a Rocky, Sandy, and 

Muddy Seabed. (Figures 32a-c). 


Figure 32a. Period 4b—Rocky Bottom 70 dB Transmission Range for a Range- 

Independent SSP and a Range-Dependent SSP. 
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Figure 32b. 


Period 4b—Sandy Bottom 70 dB Transmission Range for a Range- 
Independent SSP and a Range-Dependent SSP. 
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Period 4b, sandy bottom 70dB transmission range, range-dependent SSP (across IW front) 
(Yellow: Tx loss < 70 dB; Blue: Tx loss > 70 dB) 
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Figure 32c. Period 4b—Muddy Bottom 70 dB Transmission Range for a Range- 

Independent SSP and a Range-Dependent SSP. 
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E. PERIOD 5: IN THE WARM AND COLD ERONT OE A STEEP TIDAL 

BORE 

Period 5 was at approximately yearday 174.22, at the warm and cold interface of a 
steep-face tidal bore. Two simulation runs were conducted at data points that were 5 
minutes or 75 m apart, with one data point at the warm front of the tidal bore and the 
other at the cold front of the tidal bore. The exact time of analysis and the 70 dB 
transmission range for a range-independent SSP and a range-dependent SSP are shown in 
Figure 33 and Figure 34, with Figure 31 starting in the warm front and Figure 32 starting 
in the cold front. A sandy bottom was used in the analysis. 

While the range-dependent SSP (across internal wave fronts) did not show much 
change in the 70 dB transmission range between the two simulation time periods (8.3 km 
and 8.9 km), it was a different story for the range-independent SSP (along internal wave 
fronts). For the range-independent SSP the coming of the cold front of the internal tidal 
bore to displace the warmer waters resulted in the 70 dB transmission range increasing 

from 12.1 km to 32.1 km in just 5 minutes or just 75 m apart. 
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Figure 33. Period 5—Sandy Bottom 70 dB Transmission Range at the Warm 

Front of a Tidal Bore. 
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Figure 34. Period 5—Sandy Bottom 70 dB Transmission Range at the Cold 

Front of a Tidal Bore. 
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F. SUMMARY AND INTERPRETATION OE SIMULATION RESULTS 

The results of the transmission loss runs are summarized in Table 4. 


Table 4. Transmission Ranges for the Various Time Periods. 


Period 

Centered 

yearday 

Internal tide 

Solitons 

70 dB transmission loss range 

1 

175-177 

-Well-defined 

Lew 


Seabed 

Along IW 

Across IW 

% change 




-Small 



Rocky 

40.0 km 

21.5 km 

-46% 




temperature 



Sandy 

27.3 km 

6.2 km 

-77% 




difference 



Muddy 

26.5 km 

2.5 km 

-91% 




(2.5°C) 






2 

205-207 

-Well-defined 

Many 


Seabed 

Along IW 

Across IW 

% change 




-Large 



Rocky 

37.5 km 

18.3 km 

-51% 




temperature 



Sandy 

12.2 km 

6.5 km 

-47% 




difference 



Muddy 

1.8 km 

2.5 km 

39% 




(7°C) 






3 

210-212 

-Weak 

Many 


Seabed 

Along IW 

Across IW 

% change 




-Large 



Rocky 

36.8 km 

18.5 km 

-50% 




temperature 



Sandy 

17.2 km 

9.5 km 

-45% 




difference 



Muddy 

7.5 km 

3.7 km 

-51% 




(6°C) 






4a 

200.6- 

-Well-defined 

Lew 


Seabed 

Along IW 

Across IW 

% change 



202.6 

-Lairly large 



Rocky 

39.2 km 

18.3 km 

-53% 



Cold front 

temperature 



Sandy 

19.5 km 

6.7 km 

-66% 



oflW 

difference 



Muddy 

4.3 km 

2.0 km 

-53% 




(5°C) 






4b 

200.9- 

Same as period 

Same as 


Seabed 

Along IW 

Across IW 

% change 



202.9 

4b 

period 4b 


Rocky 

37.8 km 

19.2 km 

-49% 



Warm 




Sandy 

10.0 km 

5.8 km 

-42% 



front of 




Muddy 

2.0 km 

2.0 km 

0% 



IW 







5 

174.2- 

Same as period 

Same as 


Seabed 

Along IW 

Across IW 

% change 



176.2 
Edge of 

1 

period 1 


Sandy 

(warm) 

12.0 km 

8.3 km 

-31% 



tidal bore 




Sandy 

(cold) 

32.0 km 

8.5 km 

-73% 











Referring to Table 5, for a rocky bottom, sound traveling across internal wave 
fronts (range-dependent SSP) had its 70 dB transmission loss range reduced between 
46% to 53%, when compared to sound traveling along the internal wave front (range- 
independent SSP). 
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The along internal wave 70 dB transmission ranges were from 36.8 km to 40.0 km 
across the various time periods, giving a variability of 3.2 km. The across internal wave 
70 dB transmission ranges were from 18.3 km to 21.5 km across the various time periods, 
giving a variability of 3.2 km. 


Table 5. 70 dB Transmission Ranges for a Rocky Bottom. 


Period 

Unique characteristics 

Along IW 

Across IW 

% change 

1 

-Sound source in a well-defined bottom cold layer 
-Internal tides only 

40.0 km 

21.5 km 

-46% 

2 

-Internal tides and solitons 

37.5 km 

18.3 km 

-51% 

3 

-Solitons only 

36.8 km 

18.5 km 

-50% 

4a 

-Sound source in the cold front of an internal tide 

39.2 km 

18.3 km 

-53% 

4b 

-Sound source in the warm front of an internal tide 

37.8 km 

19.2 km 

-49% 


Variability: 

3.2 km 

3.2 km 



Referring to Table 6, for a sandy bottom, sound traveling across internal wave 
fronts (range-dependent SSP) had its 70 dB transmission loss range reduced between 
31% to 77%, when compared to sound traveling along the internal wave front (range- 
independent SSP). 

The along wave front 70 dB transmission ranges were from 10.0 km to 32.0 km, 
giving a variability of 22 km, while the across wave front transmission ranges were from 
5.8 km to 9.5 km, giving a variability of 3.7 km. It can be observed that a sound source 
placed in the cold front of an internal wave supported an extended transmission range as 
compared to a sound source placed in the warm front of an internal wave, and this 
resulted in a large variability in the along internal wave front 70 dB transmission range 


Table 6. 70 dB Transmission Ranges for a Sandy Bottom. 


Period 

Unique characteristics 

Along IW 

Across IW 

% change 

I 

-Sound source in a well-defined bottom cold layer 
-Internal tides only 

27.3 km 

6.2 km 

-77% 

2 

-Internal tides and solitons 

12.2 km 

6.5 km 

-47% 

3 

-Solitons only 

17.2 km 

9.5 km 

-45% 

4a 

-Sound source in the cold front of an internal tide 

19.5 km 

6.7 km 

-66% 

4b 

-Sound source in the warm front of an internal tide 

10.0 km 

5.8 km 

-42% 

5a 

-Sound source in the warm front of a tidal bore 

12.1 km 

8.3 km 

-31% 

5b 

-Sound source in the cold front of a tidal bore 

32.0 km 

8.5 km 

-73% 


Variability: 

22.0 km 

3.7 km 
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Referring to Table 7, for a muddy bottom, the along internal wave front 70 dB 
transmission ranges were from 2.0 km to 26.5 km (24.5 km variability), while the across 
wave front transmission ranges were from 2.0 km to 3.7 km (1.7 km variability). The 
only instance where the muddy seabed could support a longer range 70 dB transmission 
range was in period 1, which had a thin well-defined cold bottom layer. If this period was 
removed from the analysis, the muddy bottom along wave front 70 dB transmission 
ranges would be from 2.0 km to 7.5 km (5.5 km variability). 

As for percentage change in the 70 dB transmission range for the across internal 
wave fronts (range-dependent SSP) sound propagation as compared to the along internal 
wave fronts (range-independent SSP), a range gain of 39% (from 1.8 km to 2.5 km) was 
observed in period 2. This range gain only occurred when the along wave front sound 
propagation ranges were short (at approximately 2 km). The range gain could occur as 
when the grazing sound wave angle increases from beyond the Brewster angle, where all 
sound energy is absorbed by the muddy seabed, there will be some energy that will be 
reflected back into the water. All the other cases had range reduction from 0-91%. 


Table 7. 70 dB Transmission Ranges for a Muddy Bottom. 


Period 

Unique characteristics 

Along IW 

Across IW 

% change 

1 

-Sound source in a well-defined thin bottom cold layer 
-Internal tides only 

26.5 km 

2.5 km 

-91% 

2 

-Internal tides and solitons 

1.8 km 

2.5 km 

39% 

3 

-Solitons only 

7.5 km 

3.7 km 

-51% 

4a 

-Sound source in the cold front of an internal tide 

4.3 km 

2.0 km 

-53% 

4b 

-Sound source in the warm front of an internal tide 

2.0 km 

2.0 km 

0% 


Variability: 

24.5 km 
(5.5 km if 
period 1 is 
excluded) 

1.7 km 
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VIII. HOW INTERNAL WAVES AFFECT SOUND TRAVEL PATHS 


Observing that simulation runs that factored in the effects of internal waves 
generally resulted in greater transmission loss compared to simulation runs that do not 
account for the effects of internal waves, it would be logical to conclude that internal 
waves cause the sound waves to refract downward more. This leads them to strike the 
seabed at a larger grazing angle, resulting in fewer “rays” experiencing total internal 
reflection; for “rays” beyond the critical angle, it pushes them into a steeper grazing angle 
with a resulting smaller reflection coefficient. 

To verify this hypothesis, a ray trace was done for period 4a, at a launch angle of 
1° below the horizontal plane, as shown in Figure 35. The top ray trace was without the 
effects of internal waves, and the ray in the positive ranges (0 to 20 km) had 15 bottom 
bounces in the positive ranges before it exited the plot. For the lower ray trace plot, which 
accounted for internal waves, the ray in the positive ranges had 17 bottom bounces before 
it exited the plot, which is two more bottom interactions than the ray trace without 
internal waves. Comparing the internal wave ray trace plot with the SSP plot (Figure 35, 
bottom), it can be seen that when the sound wave passes through the warm front of an 
internal wave, at the 5-10 km region, with an associated larger temperature difference, 
sound waves get “depressed” more and strike the seabed at a faster rate than when 
passing through a colder mass of water. This would cause more transmission loss if the 
grazing angle is greater than the critical angle for which total internal reflection can occur 
for a rocky and sandy bottom. 

As for the ray in the negative ranges of the plot (0 to -20 km), the top ray trace 
had 15 bounces, while the bottom ray trace had 14 bottom bounces, which was one fewer 
bottom bounces. Internal waves still caused more losses in this scenario as the ray trace 
without internal waves had the ray refracted downward at a constant depth of 20 m, while 
for the ray trace with internal waves, the ray got refracted downward mostly at shallower 
depths of less than 20 m. This means that the ray that experiences internal waves would 
have had a greater grazing angle when striking the bottom, which translates to more 
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losses overall, even when it had one fewer bottom interaetions compared to the ray trace 
without the effects of internal waves. 


Figure 35. Ray Trace without Internal Waves (top), Ray Trace with Internal 
Waves (middle), and Temperature Profile (bottom). 
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In conclusion, ray tracing of sound waves confirmed the hypothesis that internal 
waves cause sound waves to refract downward more, causing more bottom interactions 
over the same distance or increasing the angle of grazing when striking the seabed. The 
net effect was more losses over the same distance as compared to a model without 
internal waves. This would explain why simulations that accounted for the effects of 
internal waves via range-dependent SSPs resulted in higher transmission loss as 
compared to simulations that did not account for the effects of internal waves. 
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IX. CONCLUSION 


The simulation of the effeets of internal waves on sound propagation, based on 
real world internal wave data, showed that internal waves introdueed great variability in 
sound propagation ranges in the shallow waters of the eontinental shelves, as the internal 
waves modulate the SSP in time seales as short as a few minutes. 

Assuming internal waves in the experiment area, which are propagating in the 
west-east direction, have plane wave structures, sound waves propagating in a 
northerly/southerly direction will then run along or in parallel to the wave fronts of these 
internal waves (see Figure 36). In this scenario, the location of the sound source with 
respect to the warm or cold front of an internal wave front has a great effect on the sound 
transmission range, with a cold front allowing a farther transmission range compared to a 
warm front. In one of the simulation runs where the warm and cold front interface of an 
internal wave was in the form of a steep tidal bore, a transmission range of 32.0 km was 
reached in the cold front of the bore, but dropped to 12.0 km in the warm front of the 
bore, which was in a data point just 5 minutes or 75 m apart. 
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Figure 36. Along Wave Front Variability in Sound Transmission 

Range Due to the Warm and Cold Fronts of Internal Waves. 
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Sound waves propagating along the warm front of an internal wave have greater 
attenuation compared to sound waves propagating along the cold front of an internal 
wave. 


For sound waves propagating in an easterly/westerly direction, they will cut 
across the wave fronts of internal waves, resulting in them having shorter transmission 
ranges as compared to sound waves of the same power level propagating in the 
northerly/southerly direction (i.e., in parallel to the wave fronts of the internal waves). 
This scenario is illustrated in Figure 37. In simulation runs conducted using a sandy 
seabed, this resulted in a 31-77% reduction in transmission ranges (from 10.0-32.0 km to 
5.8-9.5 km) when compared to sound waves that travel along the internal wave fronts. 
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Figure 37. Across Wave Front Sound Transmission 
Range Reduction Effects of Internal Wave. 
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Sound waves have greater attenuation when propagating aeross internal wave fronts 
when eompared to propagating along internal wave fronts. 


In conclusion, for sound propagation along internal wave fronts, internal waves 
cause great variability in propagation ranges, with cold fronts allowing farther 
propagation range than warm fronts. As for sound propagation across internal wave 
fronts, internal waves cause them to have a shorter propagation range compared to the 
sound waves of the same power level that are traveling along an internal wave front. 

Hence, the current practice of dropping an XBT or CDT into one spot to get the 
SSP for that location, and using that one SSP for “range of the day” estimation of 
acoustical systems performance, with the underlying assumption that the SSP does not 
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change with location, range, or time may not be the best approach in the shallow waters 
of the continental shelves. 

Going forward, in scenarios where an accurate estimation of acoustical system 
range performance is important for “mission success,” such as in the search for and 
location of an aircraft “black box” that is suspected to be on the seabed, the effects of 
internal waves on the “black box” underwater sound beacon transmission range would 
have to be accounted for in the modeling and simulation runs during “mission planning,” 
to ensure there are no inadvertent gaps in the search operation. 

This study conducted analyses of the effects of internal waves on sound 
propagation using the assumption that internal wave fronts take on the structure of plane 
waves. While this allowed a baseline generic assessment of the effects of internal waves 
on sound propagation, internal waves in the shallow waters of the continental shelves 
have been observed to have more complex structures that are closely tied to an area’s 
bathymetry (Jackson 2004). Hence, to accurately account for the effects of internal waves 
on sound propagation ranges, further data measurements are required at the locations 
where underwater acoustical operations have to be carried out. 
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